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Abstract 

The review of the mathematical treatment of plas- 
mon resonances as an eigenvalue problem for specific 
boundary integral equations is presented and general 
properties of plasmon spectrum are outlined. Promis- 
ing applications of plasmon resonances to magnetics 
are described. Interesting relation of eigenvalue treat- 
ment of plasmon resonances to the Riemann hypoth- 
esis is discussed. 



1 Introduction 

Plasmon resonances may occur in metallic (gold or 
silver) nanoparticles under the following two con- 
ditions: 1) dielectric permittivity of nanoparticles 
is negative and 2) the free-space wavelength of in- 
cident radiation is appreciably larger than geomet- 
ric dimensions of nanoparticles. These two condi- 
tions are naturally and simultaneously realized at the 
nanoscale. For this reason, plasmon resonances are 
intrinsically nanoscale phenomena. These resonances 
result in powerful nanoscale sources of light which 
have many important applications in such areas as 
nano-lithography, near-field optical microscopy, sur- 
face enhanced Raman scattering (SERS), biosensors, 
nanophotonics, optical and magnetic data storage, 
etc. 

In this paper, the eigenvalue approach [l]-[3] to 
the analysis of plasmon resonance modes in nanopar- 
ticles is reviewed and general properties of plasmon 



spectrum are presented. The excitation conditions of 
desired plasmon resonance modes are discussed along 
with the time-dynamics of excitation and dephasing 
of these modes. Some selective applications of plas- 
mon resonances to magnetics arc described. These 
applications include plasmon resonance enhancement 
of magneto-optic effects, all-optical magnetic record- 
ing and thermally assisted magnetic recording. Fi- 
nally, interesting relation of eigenvalue treatment of 
plasmon resonances to the Riemann hypothesis is dis- 
cussed. 

2 Plasmon Resonances as an Eigen- 
value Problem 

To start the discussion, we first stress that plasmon 
resonances are electrostatic in nature, because they 
occur at free-space wavelengths that are much larger 
than particle dimensions. Since the permittivity of 
nanoparticle is negative, the uniqueness theorem of 
electrostatics is not valid, and nonzero source- free so- 
lutions of electrostatics may appear. These nonzero 
source-free solutions are the manifestation of plasmon 
resonances. To find these source-free solutions (i.e., 
plasmon resonance modes), surface electric charges 
a(M) on the nanoparticle boundary can be intro- 
duced (see Figure 1). The electric field created by 
these surface charges in free space satisfies all elec- 
trostatic equations except the boundary conditions 
for the normal component of electric field on the 
nanoparticle surface S. It can be shown [1,2] that 
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Figure 1: (a) Nanoparticlc with dielectric permittiv- 
ity e; (b) Modeling of plasmon resonance modes by 
using the distribution of surface charges a(M). 



the latter boundary condition will be fulfilled if the 
density of electric charges a(M) satisfies the following 
boundary integral equation: 



°{Q) = 



2^ 



a{M) r -^^dS M , 



MQ 



where 



A 



(1) 



(2) 



e + £o 

and all other notations have their usual meaning. 

It is apparent that nonzero solution of homoge- 
neous integral equation (1) may exist only for special 
values Afe, which are the eigenvalues of the integral 
operator in (1): 



£fc + £o ' 



(3) 



As soon as eigenvalues Afc are found, the correspond- 
ing resonance values of permittivity efc can be deter- 



Figure 2: Modeling of plasmon resonances by using 
double layer charges. 

mined from (3) and then the resonance frequency ujk 
is calculated by using the nanoparticlc dispersion re- 
lation e(w): 

e fe = e'(wfc) = Re[e(w fe )] . (4) 

The eigenfunctions CTfc(M) corresponding to Xk create 
fields which can be construed as plasmon resonance 
modes. 

There is another, dual approach to the analysis of 
plasmon resonance modes. In this approach, dou- 
ble layers of electric charges r(M) are introduced on 
nanoparticle boundaries instead of single layer of elec- 
tric charges a(M) (see Figure 2). The distribution of 
these double layers must be chosen in such a way 
that they produce in free space the same field of elec- 
tric displacement D as may exist for actual plasmon 
modes in the presence of nanoparticles. It can be 
shown (see [1,2]) that the latter will be the case if 
the density r(M) of double layers satisfies the follow- 
ing boundary integral equation: 

r(Q) = ^c(r(M) rQM 3 nM dS M , (5) 



2tt 



' QM 



where A is given by formula (2). It is apparent that 
the last integral equation is adjoint to equation (1) 
and, as discussed below, they have the same set of 
eigenvalues Afe. The eigenfunctions of equations (1) 
and (5) are biorthogonal, 



j> a k {M)T t {M)dS M 



(6) 
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In numerical calculations, boundary integral equa- 
tions are discretized and reduced to eigenvalue prob- 
lems for fully populated matrices. Since these matri- 
ces are generated through the discretization of inte- 
gral operators with 1/r-type kernels, numerical com- 
putations can be considerably sped up by using the 
fast multipole method [4,5]. 

It can be shown [6,7] that the spectrum of integral 
equations (1) and (5) has the following interesting 
properties. For any shape of nanoparticles the eigen- 
values of integral equations (1) and (5) are real, which 
implies (according to the Frcdholm theory) that these 
eigenvalues are the same for both equations. More- 
over, Ao = 1 is an eigenvalue, while for all other eigen- 
values, 

|A fc | > 1. (7) 

The eigenvalue Ao = 1 does not belong to the plas- 
monic spectrum. Its respective eigenfunction <7o(M) 
corresponds to the classical Robin problem of distri- 
bution of electric charges over the surface of an ideal 
conductor. All other eigenvalues correspond to plas- 
mon resonance modes and according to (3) and (7) 
these modes may exist only for negative efe, as ex- 
pected. It can also be shown that the eigenvalues X k 
are simple poles of resolvents of integral equations (1) 
and (5), respectively. 

Unique spectral properties of plasmon resonances 
occur for metallic nanowires, i.e., in the two- 
dimensional case [1-2]. In this two-dimensional case, 
for any shape of nanowire cross-section the set of 
eigenvalues X k (k > 1) of the corresponding integral 
equation 



<t(Q) 



(8) 



consists of pairs of positive and negative eigenval- 
ues of the same absolute value. This phenomenon of 
"symmetric (twin) spectrum" is due to unique sym- 
metry properties of the mathematical formulation 
of the plasmon resonance problem which appear in 
the two-dimensional case due to the existence of the 
stream function. The twin-spectrum phenomenon re- 
sults in two distinct bands of plasmon resonances 
with relative dielectric permittivities being recipro- 
cals to one another. This twin-spectrum phenomenon 



will be used below to establish some relation of plas- 
mon resonances to the Ricmann hypothesis. 

It is apparent that the mathematical structure of 
integral equations (1) and (5) is invariant with re- 
spect to the scaling of S, i.e., scaling of nanoparti- 
cle dimensions. This implies the unique property of 
plasmon resonances: resonance frequencies depend 
on particle shape but not particle dimensions, pro- 
vided that they remain appreciably smaller than the 
free-space wavelength. 

It is clear that the integral operator in equation 
(1) is not self-adjoint in L2(S) because the kernel 
of this operator is not symmetric. For this reason, 
the eigenfunctions <7j(M) and <7fc(M) corresponding 
to different eigenvalues Aj and A& are not orthogonal 
on S in the usual sense (in i2-norm). Nevertheless, 
it can be shown that electric fields Ej and corre- 
sponding to eigenfunctions <Ji{M) and afc(M) satisfy 
the following strong orthogonality conditions: 



I 
I 



E, • E k dv = 0, 



E* • Efccfo = 0. 



(9) 



(10) 



The peculiar feature of these strong orthogonality 
conditions is that they hold separately in regions V + 
and V~. 

As previously discussed, plasmon resonances occur 
when the real part of dielectric permittivity is nega- 
tive. This makes the classical formula 



(11) 



for the time-average of stored electric energy density 
meaningless. It turns out that the last formula can be 
properly modified to be valid for slightly lossy (trans- 
parent) dispersive media. The modified formula is 



4 da 1 1 ■ 



(12) 



The presented analysis of plasmon modes is purely 
electrostatic in nature and completely ignores radia- 
tion phenomena. It is shown in [2] that the radiation 
corrections can be mathematically treated as per- 
turbations with respect to a small parameter which 
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is the ratio of particle dimension (its diameter) to 
the free-space wavelength. It is demonstrated that 
the first-order radiation corrections to the resonance 
value of dielectric permittivity (hence, to the reso- 
nance frequency) and to plasmon mode electric field 
are equal to zero for any shape of metallic nanopar- 
ticle. The explicit second-order radiation corrections 
for the resonance values of dielectric permittivity are 
derived in [2] and it is shown that in the particu- 
lar case of spherical nanoparticles these second-order 
radiation corrections coincide with the radiation cor- 
rections obtained from the classical Mie theory. 

The excitation of plasmon resonance modes by the 
incident light can be studied by using the biorthogo- 
nal expansion of electric charges <j(M,t) induced on 
S during the excitation process: 



a(M,t) = Y^a k {t)a k {M), 



k=\ 



a k (t) = <j> a{M,t)T k {M)dS M - 



(13) 



(14) 



It is clear that the time evolution of the expansion co- 
efficient a k (t) reveals the time-dynamics of excitation 
(and dephasing) of the plasmon mode corresponding 
to the eigenfunction a k (M). This dynamics has been 
studied in [3], and below, only a few, most salient, 
results from [3] are presented. Namely, the following 
analytical expression for the steady state a k ss \t) of 
the expansion coefficient in (14) can be derived in the 
case of resonance excitation: 



4 ss \t) = 

- (E • Pfe) 



e'(uj k ) - e Q 



e"(u; k ) 



cos ui k t + smojkt 



(15) 



Here, Eo is the electric field of the incident radia- 
tion, while pfc is the dipole moment of the fc-th plas- 
mon mode. The last expression clearly reveals that 
the desired plasmon modes are most efficiently ex- 
cited when the directions of their dipole moments are 
parallel to Eo and when |e'(wfc)| 3> |e"(u;fc)|, where 
e"(wfc) = Im[e(w/;)]. It is worthwhile to note that the 
ratio e'(u k )/e"(co k ) is most appreciable for gold and 
silver when the free-space wavelength is between 700 



nm and 1100 nm, and this ratio is appreciably higher 
for silver than for gold. Thus, as far as the quality of 
plasmon resonances is concerned, silver is "gold" and 
gold is "silver" . This fact has long been appreciated 
in the SERS research area where silver nanoparticles 
have been predominantly used in experiments. For 
off- resonance excitation with frequency loq =/= Lo k , the 
following expression has been derived [3] : 



a[ ss) (t) = (E • Pk)C(u ) cos(w i + 4>), (16) 



where 



C(w ) = 



[e fe -e'( Wo )] 2 + [e"M] 2 ' 



(17) 



The last formula is instrumental for evaluation of the 
sharpness (width) of plasmon resonances. 

In the case of the Drude model for the dispersion of 
dielectric permittivity, simple analytical expressions 
for time-dynamics of excitation and dephasing of par- 
ticular plasmon modes have been derived in [3] . 

3 Applications of Plasmon Reso- 
nances to Magnetics 

In this section, only a few sample applications of plas- 
mon resonances to magnetics are briefly discussed. 
We first consider the enhancement of magneto-optic 
effects in garnet films by plasmon resonances excited 
in metallic nanoparticles embedded in these films. It 
is known that on the macroscopic level magnetic gar- 
nets act as gyrotropic media that discriminate be- 
tween right-handed and left-handed polarizations of 
light. This results in the Faraday magneto-optic ef- 
fect. On the microscopic level, magneto-optic effects 
are controlled by spin-orbit coupling, which is de- 
scribed by the Hamiltonian which depends on lo- 
cal electric fields. These fields can be optically in- 
duced through plasmon resonances in metallic (gold) 
nanoparticles embedded in garnets, which may even- 
tually lead to the enhancement of magneto-optic ef- 
fects. In this way, plasmon resonances can be prac- 
tically utilized for the enhancement of the Faraday 
effect as well as for the probing of the origin of this 
effect on the very fundamental microscopic level. 
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Some experimental research [8] has been conducted 
to test the feasibility of the idea described above. 
In particular, garnet films have been grown over 
(lOO)-oriented substituted gadolinium gallium gar- 
net (SGGG) substrates partially populated with gold 
nanoparticles. These gold nanoparticles have been 
created by evaporating thin layers (about 5 nm) of 
gold on selected areas of SGGG garnet substrates 
and subsequently annealing them in air at about 830 
°C. It has been observed that these gold nanopar- 
ticles have different dimensions but mostly the same 
hemispherical shape. For this reason, they resonate 
at practically the same frequency (wavelength) due 
to scale invariance of plasmon resonances. By us- 
ing the liquid phase epitaxy (LPE) technique, thin 
magnetic garnet films of (Bi, Pr, Y, Gd)a(Fe, Gd)sOi2 
composition have been grown over gold nanoparticle- 
populated substrates. X-ray diffraction measure- 
ments have been performed and they reveal that af- 
ter LPE growth gold nanoparticles survive in crys- 
talline form. Then, Faraday rotation and magneto- 
optic loop measurements have been performed and 
meaningful (non-negligible) enhancement of Faraday 
rotation which can be attributed to plasmon reso- 
nances in embedded nanoparticles has been observed. 

Next, we shall briefly discuss the application of 
plasmon resonances to all-optical magnetic record- 
ing. It has been recently demonstrated [9] that fem- 
tosecond magnetization reversals of about 100 /im 
spots of magnetic media can be achieved by using 
only circularly polarized laser pulses. The direction 
of magnetization reversals is controlled by the helic- 
ity of circularly polarized light. It is currently per- 
ceived that these magnetization reversals occur due 
to the combination of the local laser heating of mag- 
netic media close to the Curie temperature and the 
simultaneous action of circularly polarized light as 
an effective magnetic field with the direction paral- 
lel to the light-wave vector. Furthermore, right- and 
left-handed circularly polarized light beams act as ef- 
fective magnetic fields of opposite directions. 

The all-optical magnetization switching described 
above will be technologically feasible for magnetic 
recording only if the techniques for delivery of 
nanoscale-focused circularly polarized light are de- 
veloped. The latter can be achieved by using opti- 



cally excited plasmon resonances in uniaxial metallic 
nanostructures [10-11]. The central point is that in 
the case of uniaxial metallic nanoparticles (or nano- 
holcs/apcrtures) the eigenvalues of integral equa- 
tions (1) and (5) have geometric multiplicity of two 
and for each there are two identical (up to ninety- 
degree rotation in space) plasmon eigenmodes. These 
two plasmon modes have the same resonance fre- 
quency and can be simultaneously excited by the inci- 
dent circularly polarized light. These two simultane- 
ously excited plasmon modes form a circularly polar- 
ized plasmon mode that can be instrumental for the 
dramatic nanoscale enhancement (focusing) of the in- 
cident circularly polarized light. It is worthwhile to 
point out that each of the two above-mentioned plas- 
mon eigenmodes is not rotationally invariant, while 
their combination as a circularly polarized plasmon 
mode is rotationally invariant. It is also apparent 
that the circularly polarized plasmon modes are most 
efficiently coupled to the incident circularly polarized 
light when the direction of propagation of laser light 
coincides with the symmetry axis of the nanoparti- 
cles. 

As discussed before, plasmon resonances can be 
most efficiently excited in the 700-1100 nm free-space 
wavelength range. To obtain plasmon resonances in 
this range, oblate spheroidal nanoparticles with large 
aspect ratio or nano-rings can be used. Another al- 
ternative is to use spherical or cylindrical nano-shells. 
This is because by controlling nano-shell thickness, 
plasmon resonances can be effectively shifted into the 
desired wavelength range [12]. Our numerical simu- 
lations [10-11] demonstrated that the intensity of in- 
cident circularly polarized light can be increased 750 
times in the case of silver nano-rings and 550 times 
in the case of silver nano-shells. 

Finally, we shall briefly discuss the possible appli- 
cation of plasmon resonances to thermally assisted 
magnetic recording (TAMR). This type of recording 
has been the focus of considerable technological inter- 
est and ongoing research lately. The central issue is 
the development of optical sources of nanometer res- 
olution and high light intensity. The proper profiling 
of optical spots is also an important practical issue 
because it may reduce the detrimental effect of the 
collateral heating of adjacent recorded bits. It has 
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been demonstrated [13] that plasmon resonances in 
metallic nanoparticles and perforated metallic nano- 
films hold unique promise for the development of 
such optical sources. In particular, two promising 
designs have been identified through numerous com- 
putations. They are rectangular bar-type metallic 
nanoparticles deposited on dielectric substrates and 
perforated metallic nano-films with rectangular aper- 
tures. The numerical modeling of these structures 
reveals that plasmon resonances in rectangular bar- 
type nanoparticles produce somewhat higher light in- 
tensity optical spots than the localized plasmon res- 
onances in perforated nano-films. In addition, it has 
been found that plasmon resonances in the rectangu- 
lar bar-type nanoparticles have some advantages as 
far as proper profiling of optical spots is concerned. 

4 Relation to the Riemann Hypothe- 
sis 

The Riemann hypothesis deals with the zeta function 
which for Re(s) > 1 is defined by 



(18) 



where pi are prime numbers and the right-hand side 
of (18) is the celebrated Euler product taken over all 
prime numbers. 

In [14] (its English translation can be found in [15]) 
Riemann studied the analytical continuation of the 
zeta function defined by the infinite series in (18). 
Riemann found that £(s) extends to the whole com- 
plex plane as a meromorphic function with only one 
simple pole at s = 1 with residue 1 and derived the 
functional equation for the extended C( s )- Riemann 
then introduced the entire function 

£(A) = k s -iV- s/2 r( s /2)c( s ), (19) 



1 



(i = v 7 — T), 



(20) 



where T is the Euler gamma function. The functional 
equation for ((s) implies that £(A) is an even function 
whose zeros have imaginary parts between — i/2 and 



Im(.s-) 



£4- 



M 



*Re(s) 



Figure 3: The strip < Re(s) < 1. 



i/2. Then, Riemann conjectured that all zeros of £(A) 
are real. This conjecture is the Riemann hypothesis. 

The zeta function has zeros at negative integers 
-2, -4, which correspond to the simple poles of 
the gamma function T(s/2). These zeros are referred 
to as trivial zeros. The other nontrivial zeros are 
Sfc = i + iAfe, where £(Afc) = 0. These zeros are in 
the strip < Re(s) < 1 (see Figure 3). In terms of 
Sfe, the Riemann hypothesis is stated as follows: all 
nontrivial zeros of the zeta function have real part 
equal to 1/2. 

Riemann derived the integral representation for 
£(A) and stated that £(A) "can be developed as a 
power series ... which converges very rapidly." This 
power series can be written as follows [15]: 



°° (I'm 



n=0 



where 



/oo 
H (x) (in <Jx) n dx 

3/2 ^0) 



H(x) = 4 



d_ 

dx 



dx 



1/4 



and 



(21) 

(22) 
(23) 

(24) 



n=l 
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It can be proven [15] that for x > 1 
H(x) > 

and, consequently, 

c 2 „ > 0. 



(25) 



(26) 



(a) 

-A, 



A() Aj 



A, 



1 



►A 



It is also apparent from formulas (21)-(24) that £(A) 
assumes real values for real A and the set of its zeros 



(b) 

-A, 



A. 



A, 







consists of pairs ±Afe, fc e N. Riemann also suggested — . , r . , . > 

„t +u„ <voi„ • „j „+ r„ ™ i„ iigurc 4: (a) Spectrum of integral equation (8); (b) 

spectrum of integral equation (28). 



the validity of the following product formula, 



«o)n(i-£) 



Afe 



(27) 



which was rigorously proved by Hadamard. 

It has been observed in numerical calculations (and 
then conjectured) that zeros of £(A) are simple. 

It is our intention to discuss the relation of the Rie- 
mann xi- function £(A) to the eigenvalue treatment of 
plasmon resonances in nanoparticles. Before proceed- 
ing with this discussion, it is worthwhile to remark at 
this point that the Ricmann hypothesis stated for xi- 
function defined by formulas (21)-(24) can be treated 
as a problem of complex analysis, detached from £'s 
number-theoretic origin. It is in this context that the 
following two questions can be naturally posed: 

1. If the Ricmann hypothesis is valid, then how is 
the alignment of all zeros of £(A) along the real 
axis encoded in the mathematical structure of 
formulas (22)- (24) for the expansion coefficients 

2. What are the necessary and sufficient conditions 
(in terms of C2 n ) for the alignment along the real 
axis of all zeros of an entire function defined by 
the power series (21)? 

It is attempted below to answer these questions at 
least partially. 

It has been previously discussed that A = 1 is an 
eigenvalue of the integral equation (8) while all other 
eigenvalues are real and can be grouped in pairs ±Afe, 
k G N (see Figure 4a) . These properties of spectrum 
of integral equation (8) are valid for any sufficiently 



smooth curve L. Consider the integral equation 



<r{Q) 



cr(M) 



~MQ ' 



MQ 



dl 



M, 



(28) 



where L in the integrand stands for the length of 
curve L. It can be shown that the integral equation 
(28) has the same spectrum as integral equation (8) 
except for the eigenvalue Ao = 1 (see Figure 4b) . The 
proof follows from the facts that for At ^ 1 we have 



a k (M)dl M = 0, 



while for the Robin problem (A = 1) 



(29) 



<j> a (M)dl M + 0. (30) 



Indeed, from the well known formula 



' MQ 



dl Q = 7T, (Me L), 



(31) 



follows the validity of Eq. (29) for the solutions of 
both integral equations (8) and (28) if Afe ^ 1. Thus, 
these two equations are equivalent for those Afe . How- 
ever, this equivalence does not hold for A = 1, be- 
cause equality (29) is still valid for integral equation 
(28), while nonzero solution of integral equation (8) 
exists only under the condition (30). 

Now, it is apparent that the Fredholm determinant 
-D(A) of the integral equation (28) is the entire func- 
tion that assumes real values for real A and the set 
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of its zeros consists of pairs of real numbers ±Afc, 
k G N. As discussed before, these zeros are simple 
poles of the resolvent of the integral equation (28). 
This implies that if L does not have any symmetry 
and the geometric multiplicity of all \u is 1 (i.e., no 
accidental degeneration of Afe), then A& are simple 
zeros of the Frcdholm determinant. Furthermore, it 
has been proven (see [16]) that if a kernel of an inte- 
gral equation is Holder-continuous with respect to the 
variable of integration (which is the case for equation 
(28) when L is smooth) then the Fredholm determi- 
nant is of genus zero. This implies the validity of the 
product formula for smooth L: 



D(X)=D{Q)\{ 



1 



A 

Afe 



1 



A 



(32) 



Since D{\) is an even entire function, the following 
power series can be used for its representation: 



n=0 



(-1)' 

2nl 



b 2n ^ 



In 



(33) 



It is clear from the presented discussion that for vari- 
ous curves L Frcdholm determinants of integral equa- 
tion (28) form the class of entire functions with prop- 
erties that have been conjectured or proved for the 
Riemann xi-function £(A). For this reason, the prob- 
lem can be posed to find such curve L that 



£(A) = D(X), 



which will prove the Riemann hypothesis. The last 
formula is consistent with a spectral interpretation 
of the Riemann hypothesis advanced in the Hilbert- 
Polya conjecture. It is worthwhile to remark that this 
conjecture deals with self-adjoint operator spectral 
interpretation. The integral operator in equation (28) 
is not self-adjoint in the space of square-summable 
functions because its kernel is not symmetric. How- 
ever, this integral operator has real eigenvalues and 
they are simple poles of its resolvent. These prop- 
erties are typical for self-adjoint operators. For this 
reason, it is not surprising that this operator will be 
self-adjoint in properly constructed "energy" function 
space as suggested by the orthogonality conditions 
(9)-(10). 



Indeed, consider the Hilbert space with the inner 
product 

(v,a) = <j> v(M) (J> a(P) In — !— dl^j dl M (35) 

defined on functions with zero- mean over L. It is easy 
to see that (a, a) is positive and, up to a factor, it has 
the physical meaning of energy stored in the electric 
field created by the charge distribution a over L. By 
using the relation 



ln J d_ 

L r PQ dn Q 
1 d 



In- 



dl Q = 



In— ^- — — ( In — ) dl Q , 
r MQ on Q V r PQ ) 



(36) 



which follows from the Green's formula, it can be 
shown that the integral operator in Eq. (28) is self- 
adjoint in this "energy" function space. 

From formulas (21), (33)-(34) and the known facts 
from the theory of integral equations, the following 
recurrent relations can be deduced: 

c 2n = j> B 2n -i(Q, Q)dl Q , (37) 



C2n-1 = 0, 

B (Q,M) = K(Q,M), 



(38) 
(39) 



(34) B n (Q, M) = c n K(Q, M) 



-nj^K{Q,P)B n _ 1 {P,M)dl P , (40) 
where K(Q, M) is the kernel of integral equation (28), 



K(Q,M) = 



MQ 



(41) 



It is clear that by using formulas (37)-(41) the iter- 
ated traces of the kernel of the integral equation (28) 

Q2n = <j> K 2n {Q,Q)dl Q , (q 2 „-i = 0), (42) 
K n (Q, M) = £ K(Q, P)Jf n _i(P, M)dl P (43) 
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can be sequentially found. 

If the kernel (41) can be reconstructed from its it- 
erated traces, then curve L can be reconstructed (up 
to some scaling) from this kernel. This is because 
limM-s-Q (fMQ • hqAmq) is equal to half the curva- 
ture of L at Q and a plane curve can be reconstructed 
from its curvature. 

Frcdholm determinants can be related to the iter- 
ated traces by using another relation from the inte- 
gral equation theory, 



D'(X) 



A£(A) XD(X) 



2n-2 



(44) 



which is valid for sufficiently small A. 

The last formula offers another (unrelated to the 
previous discussion) approach to the pursuit of the 
proof of the Ricmann hypothesis. Indeed, from the 
equality 



A£(A) 



X^32nA 



2n-2 



(45) 



n=l 



and from the product formula (27), we find 

?i(x^ + ^)-|«^ M - (46) 

It is clear that the Riemann hypothesis will be proven 
if it is established that all poles of the real meromor- 
phic function on the left-hand side of (46) are real. It 
is here that the Grommer theorem can be used which 
gives the necessary and sufficient conditions for the 
validity of the Riemann hypothesis in terms of expan- 
sion coefficients q n . Namely, according to the Grom- 
mer theorem [17,18], the real meromorphic function 
on the left-hand side of (46) has only real poles, if and 
only if the expansion coefficients qi n form the positive 
sequence, i.e., Hankel quadratic forms are positive, 



N 



92i+2j+272<72j > 

ij=0 



(47) 



for any N. 

To verify inequality (47), a linear functional / de- 
fined on the set of polynomials 



p(z) = J2p 2l z 2t 



(48) 



can be introduced by the formula 

m 

f(p( z )) = Yq2i+2P2i 



(49) 



i=0 



It is clear that the inequality (47) is satisfied if for 
any polynomial 



with arbitrary real ji we have 

/(7(*))>0. 



(50) 



(51) 



From formulas (21) and (45) the following recurrent 
relation can be derived: 



C2 



n+2 - Y 



k=0 



(-l) fc (2n + l)! 
(2n-2fc)! 



C2n-2fc<72fc+2- (52) 



By interpreting this relation in terms of the func- 
tional / and by using inequalities (26), we find 



C2n+2 = f(R.2n(z)) > 0, 



where 



R2n(z) = Yl 



„2fc 
'2k z ' 



(53) 
(54) 



k=0 



(2n) (-l) fc (2n+l)! 
4 = (2n _ 2fc)! ^-2k. (55) 

The set of polynomials on which the functional / is 
positive can be appreciably extended by using the fol- 
lowing property of expansion coefficients c 2 „ of £(A): 

N 

^ C2i+2j+2ma 2 i0.2j > 0, (56) 
i,3=0 

which is valid for arbitrary real a's. The last property 
follows from formulas (22) and (25). By substituting 
formula (53) into (56) and using linearity of the func- 
tional /, we find 

N 

f ( Y R 2t+ 2 j+2m- 2 (z)a 2l a2j ] > 0. (57) 
,i=o,j=o 
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It is not clear at this time if the set of polynomials 
on which functional / is positive (see inequality (57)) 
contains all "complete-square" polynomials (50). In 
other words, it is not clear if the alignment of all 
zeros of £(A) along the real axis is fully encoded in 
inequalities (25) and (56). These inequalities do not 
account (among other things) for the fast decay of 
H(x) defined by formulas (23)-(24). 

The mathematical tools used in the presented dis- 
cussion of the Riemann hypothesis are quite modest 
by modern standards. Nevertheless, it cannot be pre- 
cluded that these tools may be helpful in the pursuit 
of the elusive proof of this hypothesis. 

This research work has been supported by NSF and 
ONR. 
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